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Momentum-conserving one-dimensional models are known to exhibit anomalous Fourier's law, with a ther- 
mal conductivity varying as a power law of the system size. Here we measure, by numerical simulations, several 
cumulants of the heat flux of a one-dimensional hard particle gas. We find that the cumulants, like the conduc- 
tivity, vary as power laws of the system size. Our results also indicate that cumulants higher than the second 
follow different power laws when one compares the ring geometry at equilibrium and the linear case in contact 
with two heat baths (at equal or unequal temperatures). 
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Introduction 

Understanding the fluctuations of the flux of heat or of the current of particles through systems in their steady state is a central 
question in non equilibrium statistical mechanics. Since the discovery of the fluctuation theorem d-Sl, one knows that the 
probability distribution of these fluctuations has some symmetry properties related to time-reversal symmetry. In most cases, the 
calculation of this distribution for a given microscopic model remains, however, a challenging issue. 




FIG. 1: We measure the cumulants of the integrated current Qt of a one-dimensional hard particle gas with alternating masses 1 and mi > 1: 
fig-UJa) for an even number of particles on a ring; fig.[T]b) for an odd number of particles in contact with two heat baths (open system). In our 
simulation, we measure the flux Q t through a section S, located anywhere on the ring in fig.Q]a) an d half-way between the two heat baths in 
fig.Q]b); we a l so measure the flux Q t averaged over the whole system. 

Over the last few years, several exact expressions for the distribution of these fluctuations have nevertheless been obtained 
for diffusive systems, such as lattice gases II7I- U4I1 (for instance the one-dimensional symmetric simple exclusion process). For 
diffusive systems in their steady state, in the two geometries of fig. [T](the ring geometry in equilibrium at temperature T and the 
open system, i.e. a finite system in contact with two heat baths at temperatures T a and Tf,), the cumulants (Q") c of the flux of 
energy Q t during a long time t take, for a large system size L, the following form: 

lirn -{Q 2 t y^ ~ jA(T) ; lim -{Q 2 t n )T s * ^B n (T) for n > 2 ; lim ±(0?)?°* * jC n (T a ,T b ) . (1) 

t-VOO t L t-too t L z t-too t L 

For the ring geometry (fig. [T]a), explicit expressions of the prefactors A(T) and B„ (T) have been obtained for the symmetric 
simple exclusion process, as well as for generic diffusive systems with one conserved quantity urn- For the open case (fig. [T|b), 
the amplitudes C n (T a , T&) are also known for the one-dimensional symmetric exclusion process |7- 9] (with expressions identical 
to those which had been previously determined for disordered one-dimensional conductors IU5l flaO ; they are also known for 
generic diffusive systems with one conserved quantity 

The 1/L dependence of (Qt)° pen /t for the open case means that diffusive systems satisfy Fourier's law [20]. Together with 
the fluctuation-dissipation theorem, Fourier's law also implies the 1/L decay of (Q| ) ring /i (or of (Q%) open /t for T a = Tb). 
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It is remarkable that the L dependence (Q]i of all the cumulants is generic for diffusive systems, with a few exceptions where 
some prefactors may vanish and the decay with L is faster [10]. All these L dependencies follow from the fact that, for large 
diffusive systems, local equilibrium holds and can be treated within the macrosopic fluctuation theory (MFT), a theory of 
diffusive systems which allows one to calculate explicitly a number of properties of non-equilibrium steady states i2ll - l25ll . For 
systems with more than one conserved quantity, much less has been done so far concerning the cumulants of the current but, as 
long as the system is diffusive, one expects, from the MFT, the same L dependencies of the cumulants as inJTJ. 

Momentum-conserving systems in one dimension are known to exhibit anomalous Fourier's law Il26l - [32ll . with an average 
current varying as a non trivial power law of the system size L: 

lim l(Q t ) opcn ~ L a -^Cr{T a ,T h ) . (2) 

The exponent a is not easy to determine ll33l - i37ll . It seems to vary with the systems studied; even for a given system, numerical 
simulations or theoretical approaches do not always agree. The current consensus is that several universality classes exist 
depending on the nature of the non-linearity of the forces between the atoms I38ll39ll . 

As major numerical efforts have been already done to determine the exponent a, our goal here is to present numerical simu- 
lations, not to determine a more accurately, but rather to look at the L dependence of higher cumulants of Q t . 



I. MEASUREMENT OF THE CUMULANTS FOR THE HARD PARTICLE GAS 

The system we have decided to simulate is a one-dimensional gas I33l[35ll40ll of point particles with hard core interactions. 
It was mainly chosen for its simple dynamics: masses follow ballistic motions between successive collisions, which are elastic. 
If the particles were all of equal masses, the velocities of colliding particles would simply be exchanged, and the transport of 
energy would be the same as for an ideal gas; this is why, as in previous studies of hard particle gases, we chose here a two-mass 
system, with alternating particles of masses 1 and mi- 

For the ring geometry, we consider an even number N of masses on a circle of length L = N, initially in microcanonical 
equilibrium at fixed energy E = N and zero total momentum: the total energy and the momentum remain of course conserved 
by the dynamics. 

For the open case in contact with two heat baths at unequal (or equal) temperatures, we take an odd number N of particles in 
a one-dimensional box of size L = N, and we choose the particles closest to the boundaries to be of mass 1 . The heat baths at 
these boundaries are implemented in the following way: whenever a particle hits a boundary, it is reflected as if a thermalized 
particle was entering the system from the bath, so that the total number of particles in the system remains N. 

As the steady state is in general not known in the open case, we started our first sample with an initial condition chosen 
at random; then, the initial configuration of each new sample was taken to be the final configuration of the previous sample. 
Therefore, apart from transient effects affecting the first samples (which represent a very small fraction of the total), the initial 
configurations of our samples are typical of steady state configurations. 




2000 4000 6000 ,0 2 • 10~ 4 4 ■ 10" 4 



FIG. 2: Second cumulant of the energy flux Q t through a section and of its space average Q t on a ring of N = 800 particles (fig-Ha). When 
plotted as functions of \jt (fig.[2]b), the two ratios {Qt) /t converge, when t — > oo, to a common value ~ 0.23. 
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Through any section of our system, the flux of energy Q t is the algebraic sum of the kinetic energies of the particles crossing 
the section S during a time interval t. In the steady state, the statistical properties of Q t depend on where the section is located 
(at least in the open geometry). On the other hand, one expects (under the assumption that the internal energy cannot grow 
indefinitely) that the long time limit of the ratios (Q 7 t l ) c /t are independent of where Q t is measured. In our simulations, we 
measure for the same samples the flux of energy Q t through a fixed section S located at position L/2 in the open geometry (and 
anywhere on the ring geometry) and its integrated value Q t averaged over the whole system. 

As shown in fig. |2]a), the cumulants obtained from these two measurements behave differently at finite time, but both exhibit 
a linear growth for large t. Fig. |2]b) shows that, when the ratios (Ql) c /t are plotted versus l/t, the two sets of data converge 
to a common value in the long time limit. This was the case for all the cumulants we were able to measure: for all the results 
shown below, the procedure of fig.|2]b) was used to estimate the asymptotic values of the cumulants (Q r t l ) c /t. 

Note that the l/t convergence of fig. 2.b) can easily be understood : if the correlation function (J(ti)J(t2)) c = 
f(t2 — ti) of the energy current J(t) = dtQt decays fast enough in the steady state or at equilibrium, then (Ql) c = 
JJ Q dtidt2(J(ti)J(t2)) c — 24 J Q /(r)dr — 2 J Q r/(r)dr becomes of the form At + B in the large time limit, so that 
(QDc/t ~ A + B/t. 

Remark: all the cumulants we could measure grow linearly with time for large t. On the ring geometry, we only observed this 
growth when performing microcanonical sampling at fixed total energy E and momentum P (we took E = N and P = 0). 
When E is allowed to fluctuate (canonical ensemble) while keeping P = 0, the cumulants exhibit a faster growth {(Qf 1 ) ~ t n ). 
This is due to the conservation of E by the dynamics on the ring, which introduces non decaying current-current correlations in 
the long time limit. 



II. SIZE DEPENDENCE OF THE CUMULANTS 



Fig. |3]and|4]show the asymptotic values of (Q r t L ) c /t we obtained for 1 < n < 4. The cumulants were calculated by averaging 
over a number of samples varying from 2 • 10 8 for A = 50 to 2 • 10 6 for A = 800, for the following systems: 

• a ring of 50 < N < 800 particles with total kinetic energy A and total momentum 0; 

• an open system of 51 < A < 801 particles between two heat baths at temperatures T a = 2 and TJ, = 1; 

• the same open system between two heat baths at the same temperature T a = Tb = 1. 

In all cases, the particle density N/L was exactly 1. The mass of the heavier particles, vri2 was taken to be 2.62 (as in l33ll ): 
when we repeated some of our simulations with a different mass ratio, mi = (1 + \/5)/2 as in rt35ll . we obtained qualitatively 
similar results (not shown here). 

The fourth cumulant of Q t is of course the hardest to obtain: we were only able to measure the asymptotic value of (Qf) c /t 
accurately for N < 201 in the open case. 




Jit 10 = 



N 

FIG. 3: Asymptotic values of (Qt)c/t on a ring of 50 < TV < 800 particles of alternating masses 1 and 2.62, with total energy N and total 
momentum 0. While the second cumulant decreases as TV"" 1 , with a ~ 0.5, the fourth cumulant increases with TV. 
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For the ring geometry, while the second cumulant decays like a power law, the fourth cumulant increases with the system 
size. Hence the picture is very different from the diffusive case ([TJ: in addition to the anomalous Fourier's law, we observe an 
increase of the fourth cumulant instead of a decay. 
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FIG. 4: Asymptotic values of {Qt) c /t on an open system with 51 < N < 801 particles, with alternating masses 1 and 2.62, between 
heat baths at temperatures T a — 2 and Ti, = 1 (fig.E]a) and T a = Tj = 1 (fig.|4]b). All cumulants seem to decrease as N a ~ 1 , with 
0.25 < a < 0.4. 



For the open case, the situation looks more similar to the diffusive case : all the cumulants seem to decrease with comparable 
power laws of the system size, albeit with a different exponent than in the diffusive case (Q]). 



III. SHORT-TIME BEHAVIOR 
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FIG. 5: Short-time behavior of (Qt)c on an open system of 281 < N < 801 particles between heat baths at the same temperature T a — Tt — 
1. In a time range increasing with system size, the second cumulant behaves as on an infinite system at equilibrium (fig.|5]a); it then exhibits 
damped oscilllations on time scales of order N (fig.|5]b). 



In our simulations (fig. |2]a), we measure the whole time dependence of the cumulants (Q") c - As fig. |5]shows, the short-time 
behavior of (Q") c is independent of the system size N over a time range which increases with N. This can be interpreted as 
the fact that, over this time range, the system behaves as if it was infinite : in turn, this allows us to study the behavior of the 
cumulants on an infinite system. Fig.|5]a) shows that, for an infinite system at equilibrium, (Ql) c grows linearly with time; other 
results (not shown here) indicate that (Qf) c is also linear in time for the same system. This is another major difference with 
infinite diffusive systems ll4lll . for which all the cumulants of Q t grow asymptotically as \/t. 
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At intermediate times, all the (Q") c exhibit periodic oscillations for n > 2, with a period proportional to system size (as 
can be seen for the second cumulant on fig.|5]b): for the second cumulant, they can be fitted by an exponentially damped sine 
function J42ll . The periods can be understood from the adiabatic sound velocity c s in the hard particle gas, given by 




here, 7 = 3 (for an one-dimensional monoatomic gas), P = NkT/L = 2 (since the average energy, E/N = kT/2, is taken 
to be 1), and p = (1 + rri2)/2, so that c s = 1.82 for iri2 = 2.62. For the open system, the period is close to 2N/c s (shown on 
fig. |5]b), which is the time for a sound wave originating in S to come back to S in the same direction, having been reflected once 
on each boundary; for the ring, our data (not shown here) exhibits a period close to Njc s , the time for a sound wave starting 
from S to go around the ring. 



IV. CONCLUSION 



In this letter, we have shown that, for a hard particle gas at equilibrium and out of equilibrium, the cumulants of the flux of 
energy scale as power laws of the system size. These power laws differ from the ones expected for diffusive systems. Hence, 
our data shows that the anomalous Fourier's law of momentum-conserving systems is also characteristic of higher cumulants 
of the current. The difference between momentum-conserving and diffusive systems can also be seen for an infinite system at 
equilibrium, in which (Qt) grows as t instead of yft. 

More numerical work is certainly needed to accurately determine the exponents characteristic of the cumulants of the energy 
flux, in particular to check whether, in the open case, all the cumulants decay with the same power law; one could try to compare 
the exponents seen for the ring and for the open geometry. 

We have focused here on the fluctuations of the flux of energy, which is one of the conserved quantities of this hard particle 
gas. As momentum is also conserved, it would be interesting to study the size dependence of the fluctuations of the momentum 
flux in a similar way. 

It would also be interesting to investigate whether similar power laws of the cumulants can be observed in other momentum- 
conserving systems, such as anharmonic chains like the Fermi-Pasta-Ulam models (which are known to exhibit anomalous 
Fourier's law). The study of the cumulants could be a good test for existing theoretical approaches l43j|44ll . such as the mode- 
coupling theory [45-48]. or the Boltzmann-Langevin equation I49ll50ll . In particular, one may wonder whether, for momentum- 
conserving systems, there exists an analog of the universal ratios of the cumulants which are expected for diffusive systems on 
a ring ill- 
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